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SUMMARY 


A prediction method is found for the noise reduction through 
a cavity— backed panel* The analysis' takes into account only cavity 
inodes in one direction. The results of this analysis are used 
to find the effect of acoustic stiffness of a hacking cavity on 
the panel behavior. The resulting changes in the noise reduction 
through the panel are significant. The results of this analysis 
show good agreement with the results of the method of Guy and 
Bhattacharya (Ref. 8). The agreement with experimental results 
obtained with the test facility described in Ref. 1, however, is 
poor. 
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CHAPTER 1 


INTRODUCTION 


This report describes part of the work performed under a National 
Aeronautics nnd Space Administration (NASA) funded research project 
on the transmission of sound through general aviation airplane struc- 
tures . 

The project started on April 15, 1976, when the Flight Research 
Laboratory of the University of Kansas began work on the grant for 
NASA, Langley Research Center, entitled "A Research Program to Reduce 
Interior Noise in General Aviation Airplanes, " NASA Grant No, NSG 1301. 
The main objectives of this first phase were to develop a noise re- 
search team at the University of Kansas and to define a long-range 
follow-up research program in interior noise. 

At the end of April 1977, the ending of the preparation phase, 
the grant was extended to carry out the continuing work, which was 
defined and partially prepared during the first phase of the project. 
During this second phase a test facility was designed and built. In 
Ref. 1 the test facility is described thoroughly. The latter made it 
possible to determine the sound transmission loss characteristics 
of various structural panels and panel treatments. 

It was assumed that the receiving chamber, in the tube behind 
the test panel, was a termination which absorbed all of the sound 
passing through the test panel. In that way the receiving chamber 
would affect the transmission of sound through the test specimen in 
the same manner as an infinite space of air. However, when testing 
of panels started, several apparent anomalies were observed. One 
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problem which appeared was that the measured resonance frequencies 
of most test panels did not agree with the calculated resonance fre- 
quencies. Another point was negative transmission loss at the funda- 
mental resonance frequency. 

It appeared that the receiving chamber behind the panel behaved 
as a finite cavity and affected the panel behavior. This cavity 
effect raised the resonance frequencies and caused the negative trans- 
mission loss. In spite of the large amounts of absorbing materials 
in the tube, the receiving chamber did not have the same impedance 
as an infinite space of air. 

In this report an attempt will be made to explain the effect 
of the cavity on the panel behavior. Expressions will be derived 
for the prediction of the noise reduction through free and cavity- 
backed panels. The influence of cavity damping will also be considered. 
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CHAPTER 2 


HISTORIC REVIEW 


According to Pretlove (Ref. 2) there are two types of panel cavity 
systems. The first type is a system in which the effect of the cavity 
on the panel behavior is negligible. The panel stiffness is consid- 
erably greater than the acoustic stiffness, and the fundamental reso- 
nance frequencies are not affected, by the cavity behind the plate. 

The second type is a system in which the effect of the cavity on the 
panel is considerable. This happens when the acoustic stiffness of 
the cavity is equal to or greater than the stiffness of the panel. 

The panel mode shapes and natural frequencies change considerably. 

For this system it is necessary to include the effect of the cavity 
in the analysis of the noise transmission through the panel in the 
cavity. 

The manner in which the acoustic cavity modes affect the flexible 
panel depends very much on the length of the cavity. When the length 
of the cavity is greater than the larger panel dimension, two acoustic 
cavity depth modes will occur before any other. 

The first acoustic cavity mode is called the "open-ended" mode. 

The depth of the cavity is one quarter of the acoustic wavelength. 

In this case, the cavity affects the panel as if there were no cavity 
at all. The acoustic stiffness due to the cavity has dropped to 
zero at this frequency. The second acoustic cavity mode, which 
occurs at about twice the frequency of the "open-ended" mode, is 
the "closed-end" acoustic depth mode. In this case, the acoustic 
stiffness approximates minus infinity. 
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When the in vacuo value of the panel natural frequency (the 
natural frequency of the panel when hacked by an infinite space) is 
less than the "open-ended” acoustic mode frequency, then the cavity 
acts as a stiffness and the natural frequency of the panel is increased- 

The second possibility is that the in vacuo panel resonance fre- 
quency falls between the "open-ended" and the "closed-end" acoustic 
modes. Then the cavity acts as a negative stiffness, and the panel 
fundamental frequency is decreased. 

Xf the depth of the cavity is less than one— half of the larger 
panel dimension, then the effect of the transverse acoustic modes of 
the cavity becomes more important. 

If the depth of the cavity becomes less than one-quarter of 
the larger panel dimension, then the lowest acoustic mode in the 
frequency scale' is the transverse acoustic ►mode and infinite values 
of acoustic stiffness occur before zero values. 

The problem of a cavity-backed panel has already been studied 
by several persons. One of the first articles about this subject 
was written by Dowell and Voss in 1963 (Ref. 3). It treats 
a case in which the cavity acts as a stiffness to the fundamental 
panel mode. The effect of the cavity on the second (anti- symmetric) 
panel mode is that of a negative stiffness. According to this theory 
the cavity effect becomes more noticeable with decreasing cavity 
length. 

Lyon (Ref. 4) treats the noise transmission problem for 
the situation in which the panel vibration is not affected by the 
cavity. The acoustic stiffness because of the cavity is assumed 
negligible, compared with the stiffness of the panel. The usual 
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practical situation is considered: the in vacuo panel fundamental 

frequency is lower than the first acoustic cavity mode. The noise 
reduction study is divided into three frequency ranges: at "low" 

frequencies, where both the panel and the cavity are stiffness con- 
trolled; at "intermediate" frequencies, where the panel is in the 
resonance controlled region and the cavity is still stiffness con- 
trolled; and at "high" frequencies, where both the panel and the 
cavity are in the resonance controlled region. 

Pretlove discusses, in Ref. 2 and 5, the free -and .the forced 
vibrations of a panel-backed cavity. The model used by Pretlove 
considers a cavity with all walls perfect acoustic reflectors, 
except for the flexible panel. The solution derived in Ref. 2 is 
exact. 

The problem of a room- or cavity-backed panel has also been 
•studied by Kihlman (Ref. 6) and Bhattacharya and Crocker (Ref. 7) 
using models somewhat similar to the ones used by Pretlove (Ref. 2 
and 5) and Dowell and Voss (Ref. 3). 

Guy and Bhattacharya (Ref. 8) discuss the influence of a finite 
cavity on the transmission of sound through a panel backed by that 
cavity. The model considered consists of an acoustically hardwalled 
rectangular cavity having at one face a flexible vibrating panel. 

The same model is also considered by the above-mentioned authors. 

The effect of panel damping can be introduced. Cavity damping, 
however, is not considered in this model. The theoretical predictions 
show, good agreement with experimental results and with Pretlove’ s 
results (Ref. 5). 
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One .of the latest articles about sound transmission through 

cavity-backed panels is Ref. 9. Dowell has developed a theoretical 

model for interior sound fields- which are created by flexible wall 

motion resulting from exterior sound fields. Included in the model 

are- the mass, stiffness and damping of the panel, as well as the 

mass, stiffness and damping (due to absorbing walls inside the cavity) 

> 

of the acoustic L cavity. Comparisons of the theoretical predictions 

with experimental results show good agreement. Contrary to Guy and 

• 4 - < ‘ 

Bhattacharya (Ref. 8), however, Dowell does not compare experimental 
noise reduction data with theoretical predictions. 

Barton (Ref. 10) presents some experimental findings on the 
effects of panel stiffness and receiving space absorption on low- 
frequency sound transmission through panels into closed spaces. Also 
a simple method for predicting the low— frequency noise reduction of 
a panel backed by a closed, absorbent cavity is presented. The 
analytical model used by Barton is analogous to that described in 
Ref. 11 (page 82), except for the absorbing backing wall. The model 
uses only one panel mode and consists of a rigid, infinite panel 
supported on springs and backed by a wall having a frequency-dependent 
impedance. The theoretical predictions show fair agreement with 
the experimental data. 

In Ref. 12 an attempt is made to predict the effect of an ab- 
sorbing cavity on the panel dynamics. Although the results show 
good agreement with experimental results, a* number of comments must 
be made. The model used is the same as the one derived in Ref. 8. 

To account for the absorbing materials in the receiving chamber, a 
resistive part is added to the cavity impedance, which in Ref. 8 
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only consists of a reactive part. However, the absorbing materials 
also cause a change in boundary conditions, and the equations as derived 
in Ref. 8 (and used in Ref. 12) do not satisfy those boundary conditions. 
In Ref. 12 an effective cavity length has been used which decreases 
with increasing- frequency. Physically this is not quite right because 
with increasing frequency the absorbing materials become more and more 
effective, and less and less sound will be reflected by the walls. This 
means that the cavity becomes larger with increasing frequency; and at 
high frequencies no sound at all will be reflected by the walls, and 
the cavity behaves like an infinite space. Hence, the cavity must 
become larger with increasing frequency instead of smaller as suggested 
in Ref. 12. 

In the following chapter, the transmission of sound through a free 
panel will be discussed. 
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CHAPTER 3 


NOISE REDUCTION THROUGH A FREE PANEL 

In this chapter the transmission of sound through a panel which 
is backed by an infinite space will be discussed. This chapter will 
be an introduction to the one in which the effect of a finite space 
(cavity) behind the panel will be discussed. 

3.1 Transmission of Sound through Panel 

When a progressive plane wave in a fluid medium impinges on the 
boundary of a contiguous second medium, a reflected wave is generated 
•in the first medium and a transmitted wave in the second medium. Here 
the boundary between the two mediums consists of a thin, flexible panel. 
The model used is sketched in Figure 3.1. 


x 



8 


The panel is assumed to be very thin. The length of the sound wave, 
is much greater than the panel thickness, h. In that case any 
wave transmission through the panel can be ignored. An acoustic imped- 
ance at z - 0 is the only effect of the plate on the acoustic waves. 

The panel at z ~ 0 forms the boundary between medium I, of charac- 
teristic impedance c^, and medium II of characteristic impedance 
^2 c 2’ w ^ ere P represents the undisturbed equilibrium density of the 
medium and c the speed of sound of the medium. Consider an incident 
plane wave traveling in medium X in the positive z-direction. This 
incident wave may be represented by: 


— . i (a). t - k. z) 

p,=Ae J r i 


(3.1) 


where A is a real constant representing the pressure amplitude of 
the wave, uk is the angular frequency and: 


k .-=! 

1 C„ 


(3.2) 


Upon striking the panel which is for convenience at z 
reflected wave: 


= 0, a 


— TT 1 (<l> t + k z) 
p = Be J r r 
r r 


(3.3) 


and a transmitted wave: 


- C a ~ k t z) 


(3.4) 


are produced. The reflected wave has a frequency a complex pressure 
amplitude B and a wavelength constant: 


to 

k - — 
r c-, 


(3.5) 
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( 3 . 6 ) 


The transmitted wave has a frequency of a complex amplitude C 
and a wavelength constant: 



There are two boundary conditions that must be satisfied at all 
times and at all points on the panel surface* The conditions are 
continuity of pressure and continuity of particle velocity at z = 0. 
These two conditions can be combined and become instead a condition 
of continuity of their ratio: 


u. + u 
r r 



( 3 * 7 ) 


at z =0, where z represents the normal specitic acoustic panel 
impedance* The incident, reflected and transmitted wave are assumed 
to be plane progressive waves. The sum of the incident and reflected 
wave can be a standing wave, however. For plane progressive waves 
the following expressions are valid for the particle velocity of the 
incident, reflected and transmitted wave, respectively: 



u 

r 



p i c i 


u 


t 



P 2 C '2 


( 3 . 8 ) 

( 3 . 9 ) 

( 3 . 10 ) 


The normal specific impedance 


z 

P 


of the panel may be written as 


follows: 
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(3.11) 


z 

p 


r 

P 


+ j ^ 


P 


where r represents the panel resistance and x the panel reactance. 

P P 

The normal specific acoustic impedance is defined as the complex ratio 
‘of the driving acoustic pressure to the resulting velocity of the 
plate surface. Substitution of the Eqs. 3.8, 3.9, 3.10 and 3.11 in 
Eq. 3.7 results in the following expression: 


p . + p 

*r 


- - p i c i 

Pi - p r 


r p + 3 x p + p 2 c 2 


(3.12) 


or: 



(r p + p 2 C 2 ~ P 1 C l ) + j X P 
(r p + p 2 c 2 + P;L c x ) + j x p 


(3.13) 


By taking the magnitude of both sides of Eq. 3.13 the following 
expression can be derived: 


B 

A 


(r p + p 2 c 2 ~ P 1 C l )2 + V 

(r p + p 2 c 2 + P 1 C l )2 + X p 2 


(3.14) 


See Appendix A for the complete derivation of Eq. 3.14. 

Substitution of Eqs. 3.1 and 3.3 in Eq. 3.13 results in 
(at z = 0) : 


B e J + V 
A a J “i c 


(3.16) 


r p 2 + 2r p p 2 c 9 + (p 2 c 2 ) 2 - ( Pl c 1 ) 2 +x t ^ 2 + j 2 x^ p 7 c„ 


'1 ''l' 


K. Pi 
P 1 1 


(r p + p 2 C 2 + P 1 C l )2 + X p 2 
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When the assumption is made that the frequency of the incident wave 
is identical to the frequency of the reflected wave: 


0) . 
x 


= (jo. 


then the phase angle 0^ can be derived from Eq- 3.16 and may be 
written as: 


0 = 
r 


arctan 


2 X P p i c i 


( r p + P 


2 c 2 ) 


- (p x c x ) 2 +x p 2 


(3.17) 


At z = 0 also the following boundary condition must be satisfied: 
the particle velocities normal to the surface are equal. This is true 
when the panel damping is zero. This boundary condition can be 
written as follows: 

u. + u = u (3.18) 

x r t 


Substitution of Eqs. 3.8, 3.9 and 3.10 in Eq. 3.18 results in: 




r 


P 2 C 2 


(3.19) 


When the assumption is made: 


to. = ( 1 ) = 0 ) = 0 ) 

x r t 


Eq . 3.19 becomes : 


A - B = 


P 1 C 1 F 
P 2 C 2 


(3.20) 
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Transmission loss is defined as 10 times the logarithm, to base 10', 
of the ratio of sound power incident on the panel to the sound power 
transmitted by the panel (Ref. 13). Transmission loss is a useful 
analytic concept, but its measurement is hampered by the lack of 
acoustic watt meters. Sound power can be derived from measurements 
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of sound pressure only under certain limited circumstances, e.g. when 
the acoustic waves are plane. The formation of a reflected wave with 
the same frequency as the incident wave will generate a pattern of 
standing waves in the fluid. These standing waves make it difficult 
if not impossible to measure the pressure amplitude of the initially 
incident plane progressive wave. 

Instead of transmission loss, noise reduction is .a more useful 
term. Noise reduction is defined as the difference in sound pressure 
level at the source side and the receiver side of the panel (Ref. 13).: 


NR = SPI^ - SPL 2 (3.26) 

where SPI^ and SPI^ represent the sound pressure level at the source 
side and the receiver side of the panel, respectively. The sound 
pressure level SPL is defined as: 


SPL = 10 



(3.27) 


where P is the measured effective pressure of the sound wave (root- 
e 

mean-square pressure) and P^^ is the effective reference pressure. 

-5 2 

A reference pressure ? ref a 2 x 10 N/m is commonly used for com- 
puting sound pressure levels in air. The effective pressure amplitude 
may be written as: 


P 

e 


P 

fi 


(3.28) 


where P is the peak pressure amplitude. Now Eq. 3.26 can be written 
as follows : 
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NR = 10 log 


(3.29) 



where P^ and P£ are root-mean-square (rms) pressure amplitudes at the 
source side and the receiver side, respectively. 

At the source side at a fixed position z^, the acoustic pressure 
p^ is the sura of the real part of the incident wave and the reflected 
wave: 


p^ = A cos (oKt - k^z^) + B cos (co^t + ^ r z i + 0 r ^ (3.30) 


where 0^ is a phase angle defined as: 


j 0 r 

B = B e 


(3.31) 


This angle 0^ measures the amount by which the reflected pressure 
leads or lags the incident pressure. To obtain the rms pressure, 
first the square of p^ (Eq. 3.30) has to be taken: 

2 A 2 + B 2 . k * C0S 2 C V ' hV . 
h — + 2 + 

9 

B cos 2 (oj t + k +0 ) 

£ Li £_ + (3.32) 

2 

2 A B cos (ok t - k^z^) cos (aj^t + k^z^-f ® r ) 


2 

Secondly, the time average of p^ (Eq. 3.32) has to be taken. The 
time average of a function is defined as: 


< f > 


lira, 
T + “ 


+ T 



f (t) dt 


(3.33) 


15 



The time average of each of the single cosines in Eq. 3.32 is zero. 

It can be shown that the time average of the cosine product is also 
zero when the frequencies ok and are not identical (see Appendix B). 

The rras pressure for a combination of two tones of different 
frequencies is : 


P 


1 



(w i * “r> 


(3.34) 


Under the condition <a. = to = u. the rms pressure amplitude becomes: 

i r I 



A 2 + B 2 + 2 A B cos (6 + 2 k- z . . 

r 1 l; 


(3.35) 


This is shown in Appendix C. 

At the receiver side (position Z 2 ) of the free panel, only the 
transmitted wave exists. The acoustic pressure p 2 is the real part 
of the transmitted wave: 


P 2 = C cos (to t t - k t z 2 + 9 fc ) (3.36) 

where 0 represents the phase difference between the incident wave 
and the transmitted wave or: 


C = 



(3.37) 


The rms pressure amplitude at the receiver side can be derived and 
is: 


P 


2 


_C_ 

■J2 


(3.38) 
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Substitution of Eqs'. 3 *35 and 3.38 in Eq. 3.29 results in the following 
expression: 


NR = 10 log 


+ B 2 + 2 A B cos (0 + 2 k- z.) 

r 11 


(3.39) 


In the case ok = = to^. The expression may also be written in the 

following manner: 


NR = 10 log 


— R R 

1 + — Tp + 2 — cos (8 r + 2 z^) 

A 


(3.40) 


where the ratio B/A follows from Eq. 3.14, C/A follows from Eq. 3.23 
and the phase angle 0^ from Eq. 3.17. In Eq. 3.40 z^ indicates the 
source microphone position. The noise reduction of a free panel is 
a function, of the panel resistance, panel reactance, frequency of 
the sound, the source microphone position and the density and speed 
of sound of the air at both sides of the panel. 
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3.2 Panel Impedance of a Free Flexible Panel with Simply Supported 
Boundary Conditions 

The differential equation of forced motion of an isotropic flat 
plate is obtained from Ref. 14 and may be written as: 

D 7 2 V 2 w y , t ) + m 92 - w (x I y> = P, (x, y, t) (3.41) 

s at 2 z 

where x and y are Cartesian coordinates in the plane as indicated by 
Figure 3 .l* The complex flexural rigidity D can be computed with the 
following expression: 

— 3 3 

. o = Ul = (1 + jn) - D (1 + jn) (3.42) 

12(1 - v 2 ) 12(1 - v 2 ) 

where n = the internal loss factor of the plate 

h = plate thickness [m] 

v = Poisson's ratio of plate material 

2 

E = Young's modulus [N/m ] 

2 

= mass of plate per unit area [kg/m ] 
v = plate deflection [m] 

2 

== forcing pressure [N/m ] 

The first term on the left-hand side of Eq. 3*41 represents the 
structural stiffness of the plate. The second term at the same side 
represents the structural inertia. On the right-hand side of Eq. 3.41 
there is the pressure loading. 

Aluminum, steel and plexiglass have, according to Ref. 13, a 
very low value for the internal loss factor q. The internal loss 
factor for this type of materials is lower than 0.01. Here the 
assumption will be made that the effect of the structural damping 
is negligible. 
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3.2.-L rree viDracion or axmpry auppuiceu Flate 

For the case of a freely vibrating plate, the external force, 
P z> is zero and the differential equation of motion becomes: 


D V 2 V 2 w (x, y, t) + m - — W ~ ^ = 0 


(3.43) 


3t 


where V 2 V 2 represents the two-dimensional Laplacian operator: 


V 2 7 2 = -2-j- + 2 
3x 4 


2 2 4 
3x 3y 3y 


(3.44) 


The solution of Eq. 3.43 can be put in the following form: 

w (x, y, t) = W (x, y) 3 (t) (3.45) 

W (x, y) = X (x) Y (y) (3.46) 

The time-dependency of the displacement 3 (t) is assumed to be 
harmonic : 

3 (t) = e jaJt (3-47) 

The solution w (x, y, t) must satisfy the boundary conditions of 
the motion at t = 0. 

These initial conditions of the plate are: (w) ^ ^ and (w) ^ _ q 

Substitute Eq. 3.45 into the differential equation of motion Eq. 3.43 
and use primes to denote differentiation with respect to the variables- 
x and y, while dots indicate differentiation with respect to time t. 
Now the governing differential equation of free vibration becomes: 

D X ,,f, (x) Y(y) $(t) + 2 D X' '(x) Y f ' (y) 6(t) + 

D X(x) Y’ ' ' ’ (y) S(t) - m g to 2 X(x) Y(y) 3(t) = 0 (3.48) 
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m X Y 


X""T + 2X"Y" +XY' 


(3.49) 


The boundary conditions are assumed to be simply supported; both the 
displacement and the edge moment are zero: 


(w) n =0 

x = 0, x = a 


3x 2 


x = 0, x = a 


(3.50) 


(w) n . = 0 ( — ) 

w 7 = 0, y - b *2 


y = 0, y = b 


The solution of the governing differential equation Eq. 3.43 is ob- 
tained by Navier’ s method as follows: 


W(x, y) = X(x)Y(y) = 2 2 W sin HH2E. s ^ n (3.51) 

a=l n-1 m a b 


which expression satisfies all the above stated boundary conditions. 
Now Eq*. 3.49 may be written as follows: 


44 2222 44 

, m tt 0 m u n u . n it ^ 

1 4 + Z 2 2^4 * 

a abb 


m s 2 

T w “ 0 


4 m 2 n 2 . m s 2 
ir (— + -2> --5T“ 

a b 


which leads to: 


r 2 / =£ * si JL 


01. = » IF 2 ( — + — ) U f- 

*2 K 2 V 


where m = 1, 2, 3, • * * and n = 1, 2, 3, 
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From Eq. 3.54 the- natural frequencies of a simply supported plate can 
be obtained. 

3.2.2 Forced Vibration of Simply Supported Plate 

The governing differential equation of motion is: 

D V 2 V 2 w (x, y, t) + m — ■ — y> — = p (x, y, t) (3.41) 

5 at 2 

The assumption will be made that the forcing pressure on the plate 
is not a function of x and y, but only a function of time* The 
forcing pressure is the sum of the incident, the reflected and the 
transmitted sound pressure: 

P z (t) = ?. + p r - ? t . (3.55) 

The differential equation of motion is linear* This means that 
the differential equation of motion can be solved separately for p^, 
p^ and p t * Except for this division, the solution of Eq. 3.41 can 
also be divided into two parts: 

w (x, y, t) = w H (x, y, t) + w p (x, y, t) (3.56) 

where w.. represents the solution of the homogeneous form of the equation 
n 

(Eq. 3.43) while Wp is the particular solution of Eq. 3.41. The homo- 
geneous solution is associated with the free vibration of the plate at 
its natural frequency. Here it is assumed that the free vibration is 
successfully damped; consequently, only the particular solution is 
considered. The total steady-state deflection of the plate is: 

w_ == w. + w + w„ (3.57) 

P i r t 
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wnere w^, and represent the complex deflection caused by the 
incident wave, reflected wave and the transmitted wave, respectively. 

In accordance with Navier ! s solution method, the particular 
solution of Eq. 3.41 can be written as follows: 


w p (x, y, t) 


e 


j0)f 2 
m=l 


co 

- T , * isnrx . mry 

S W sxn sxn -- 

_ mn a b 

n-1 


(3.58) 


which satisfies the boundary conditions of Eq. 3.50. The lateral 

load p (x, y, t) can also be expanded into a double sine series: 
2 


p (x, y, t) 

Z 


= e 


Ja,t I E 
m-1 n==l 


P 

mn 


sin 


IIHTX 

a 


sin 


mry 

b 


C3.59) 


Substitution of Eqs. 3.58 and 3.59 in the governing differential 
equation of motion, Eq. 3.41, yields: 


D W 
ran 


4 4 
m it 


+ 2 


2 2 
m ir 


2 2 
n ir 


4 4 
n ir 


] -m 


2 


0) 


W = P 
mn mn 


(3.60) 


This expression can also be written as follows: 


W = 
mn 




(3.61) 


Substitute Eq. 3.54 in Eq. 3.61 and the result is: 


W 


P ' 
mn 


mn 


, 2 2 * 
m (to - a ) 
s mn 


(3.62) 


According to Ref. 14 the coefficient of the double Fourier expansion 
of the lateral load is : 
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a b 

P = — -t* / / P (x, y) sin sin d x d y (3.63) 

mnab'' ,r z ' a b 

0 0 

Here tbe lateral load is assumed to be not a function of x and y but 
a constant: 

P z (x, y) = p Q (3.64) 

From Eqs. 3.63 and 3.64 the following result can be derived: 


16 P o 

— both m and n are odd 

P = < tt“ ran 
ran 

, 0 otherwise (3.65) 

As mentioned above, the lateral load can be divided into three separate 
loads: the incident, the reflected and the transmitted pressure wave. 

At z = 0: 

p (t) = p (3.66) 

z r o 

or: 

jw.t _ ju> t _ jrn t 

p (t) - A e + B e -Ce C3.67) 


which expression may also be written in the following manner: 

jw t j (to t + 0 ) 3 (w t + 0 ) 

P z (t) - A e + B e - C e (3.68) 

where 9^ and 0^ are the phase angles between the incident and the 
reflected wave and the incident and the transmitted wave, respectively ♦ 
In the case 'that the incident wave excites the panel as only 

one: 
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From this 


. „ jo).t 

30 )t . x 

p e = A p 
o 


expression follows: 


(3.69) 


P 


o 


- A 


(3*70) 


( 0 = 0 ). 

x 

Secondly, the reflected wave is the only wave which strikes the panel; 
then: 


. . 3 (o) t + 0 ) 

iwt ^ J r r y 
p e J = B e 
o 


or: 


p cos cot = B cos (o) t + 6 ) 
r o r r 


p sin cot = B sin (to t 4* 0 ) 
*0 r r' 


From Eqs. 3.72 and 3.73 follows: 


(0 = 0) +0 /t 

r r 


V B 


(3.71) 


(3.72) 

(3.73) 


(3.74) 


Third and last case is that the transmitted wave is the only wave 
which hits the panel: 


J»t = _ 


C e 


(t^t + 


V 


(3.75) 


or: 


o> = co^_ + 0^/t 
P o = -c 


(3.76) 
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Substitution of Eqs. 3.58, 3.70, 3.74 and 3.76 in Eq. 3.57 


results in: 


J x „ ^ TT . mux . nuy . 

w = e Z Z W sxn sxn * + 

P - t mn. a b 

m=l n=l x 


i (w t + 9 ) 03 «° 

J r r „ . mux . nuy t 

e Z Z W sin sm -r*- + 

. - mn a b 

m=l n=l r 


where : 


)<v + V 


I E W sinBS-sin^ 

a=l n=l m t a b 


(3.77) 


mn. 2 ,2 2. 

i m ir m n uo -os.; 
s mn x 


(3.78) 


mn 2.2 2. 

r mirmnCco - w ) 
s mn r 


(3.79) 


- 16 C 

mn. 2 ( 2 2. 

t m it m n (w — ca^_ ) 
s mn t 


(3.80) 


which expressions follow from Eqs. 3.62, 3.65, 3.70, 3.74 and 3.76. 


3.2.3 Normal Specific Acoustic Impedance of Panel 

The normal specific acoustic impedance z^ of a plate is defined 
as the ratio of the driving acoustic pressure to the resulting velocity 
of the plate: 


p. + p - p. 
*x r r r t 


(3.81) 


25 



where represents the panel velocity.. The panel velocity can be. 
derived from Eq. 3.77 and may be written as follows: 


Wp = j w. + j u r w r + j w t w t 


(3.82) 


Now the panel impedance can be written: 


j£«Kt jOop + ep + 6 t ) 

I = - tJL g — Z-S-S. C3.83) 

^ i (to. w. + to w + ox w. ) 

J x x r r t t 


Substitution of Eqs- 3.77, 3.78, 3.79 and 3.80 in Eq. 3.83 results in 
the following expression: 


jco.t j(u tr0 ) jCut + 0) 

; 1 A e 1 + B e - C e C . 

i 16 . " " A j(<o t + « ) « » 3(0) +6 ) » » „ 

■ ' iJ: V {<xe 2 2 P — s-i). +(ae r 2 Z r-5 — =-$ -we 2 2 — -} 

_ 2 i , , 2 2. mn r - - , 2 2 S Y mn t - - , 2 2v ran 

n m-1 n-1 (u -to. ) m=l n=l (a> - w ) m=l n=*l(o> - w ) 

s mn i nm r tnn t 


(3.84) 


where : 


, 1 . m7rx . niry 

ib * — sxn sin — 

mn mn a b 


(3.85) 


To obtain a useful expression, the following assumption will be made: 


ok = d) r = £o t = (0 


In that case Eq. 3.84 simplifies to the following equation: 


(3.86) 


i 16 „ „ sin (mix/ a) sin (mry/b) 

2 co Zj T. 2 ^ 

m ir, m=l n=l m n (to — to ) 
s mn 
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In Section 3.1 the panel impedance is defined as: 


z 

P 


r 

P 


+ jx 


(3.11) 


Combination of Eq. 3.11 and 3.86 leads to the following result: 


r = 0 
P 


x = 
P 


- m ir^/16 oj 
s 


sin (mnx/a) sin (rary/b) 


2 2 ' 
m=l n=l 


' 2 2 
in n (a) - (d ) 

mn 


( 3 . 88 ) 


The panel resistance r^ is zero. This results from the assumption 
made earlier, that the internal damping of the plate is zero. 
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CHAPTER 4 


NOISE REDUCTION THROUGH A CAVITY-BACKED PANEL 


In the preceding chapter the transmission of sound through a 
panel backed by an infinite space has been discussed. In this chapter 
the effect of a finite space behind the panel will be discussed. As 
will be shown, this effect is quite noticeable, especially in the case 
of a panel with a low stiffness.. 


4.1 Transmission of Sound through Panel 


The model which will be considered in this chapter is sketched in 
Figure 4.1. 



Figure 4.1: Sketch of Panel and Cavity Geometry 


At z = 0 the side of the box consists of a flexible panel which has 


a normal specific acoustic impedance z . At z = H the rear wall of 






the cavity is situated, and this wall has a normal specific impedance 
z . The model used only considers sound waves in one direction: the 

z-directions . Thus cavity modes in the x-direction and the y-direction 
will he omitted. This assumption is reasonable as long as the length 
(depth) of the cavity is much larger than the width and height of the 
cavity. 

Again there is an incident sound wave in medium I which may be 
represented by the following expression: 


_ j (w . t - k z) 

p. = A e 

X 


( 4 . 1 ) 


where A is the amplitude of the sound wave and the wave length constant 
is : 


k. 

x 



( 4 . 2 ) 


When the sound wave strikes the panel, a reflected wave and a trans- 
mitted wave will be produced. The reflected wave may be written as 
follows: 


3 (a) t + k z) 
— — r r 

p = B e 
*r 


( 4 . 3 ) 


where B represents the complex amplitude 
and the wavelength constant is: 

03 


k 


r 


r 


c 


1 


of the reflected sound wave 


( 4 . 4 ) 


The transmitted wave is: 

_ _ jkt - k z) 

P t = c e 


( 4 . 5 ) 
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where C is the complex amplitude of the wave and: 


k 


t 



(4.6) 


Now the panel is not backed by an infinite space, as in Chapter 3. 
Because of the backing wall, a reflected wave will be produced when 
the transmitted wave strikes the backing wall. This transmitted- 
reflected wave may be represented by the following expression: 



— j ( “tr c 
D e 


k z) 
tr J 


(4.7) 


where : 



(4.8) 


and D is the complex pressure amplitude of the transmitted-reflected 
wave . 


The normal specific acoustic impedance of the wall at z = & is: 


z 

w 


r + j 
w J 


w 


(4.9) 


where represents the resistive component and x^ the reactive com- 
ponent of the wall impedance. Here it will be assumed that the wall 
is very stiff and does not vibrate, in which case x - 0. When the 
wall is a perfect reflector, the impedance of the wall z - 
However, when the wall also consists of absorptive materials, then 

z = r # 00 ♦ In this section the wall impedance will be written as: 
w w » 


z 


w 


w 


(4.10) 
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The separate boundary conditions of continuity of pressure and 
continuity of particle velocity at z — £ become a condition of conti- 
nuity of their ratio: 


+ ?tr - 
= z 


u t + u tr 


w 


(4*11) 


at z = h In Eq, 4*11 and u fcr represent the particle velocity 
of the transmitted wave and the transmit ted- reflected wave* respectively'. 
The transmitted wave and the transmitted-reflected wave both are assumed 
to be plane progressive waves . In that case the particle velocity of 
the transmitted wave may be written as: 


U t p 2 c 2 


(4.12) 


while the particle velocity of the transmitted-reflected wave is: 


— tr 

U 

tr P 2 C 2 


C4.13) 


Substitution of the Eqs. 4.10, 4.12 and 4.13 in Eq. 4.11 results in 
:he following expression: 


P 4- p 

H tr 


^t p tr 


w 


P 2 c 2 


(4.14) 


?rom Eq. 4.14 the following expression may be derived: 


tr 


r w - p 2 C 2 
r w + P 2 C 2 


(4.15) 
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Substitution of Eqs. 4.5 and 4.7 in Eq. 4.15 leads to the following 
result : 


D e 


j (w. t + k Z + 9 ) 

J tr tr tr' 


j (u._t - k.& + © ) 
J t t t ' 


C e 


r w ~ P 2 C 2 
r w + P 2 C 2 


(4.16) 


where 0^^ is the phase angle between the transmitted-reflected wave 
and the incident wave and 0^ represents the phase angle between the 
transmitted wave and the incident wave. By taking the magnitude of 
both sides of Eq. 4.16, the following result is obtained: 


D _ r w P 2 C 2 
c " r w + P 2 c 2 


(4.17) 


According to Eq. 4.15 the ratio p / p fc is real and does not 
have an imaginary part. Eq. 4.16 may be written as follows: 


D (cos a + ,j sin a). _ r w P 2 C 2 
C (cos 8 + j sin 6) r w + p 2 c 2 


(4.18) 


where : 


a = 


io. t + k. Z 
•tr tr 


+ 8 


tr 


6 =» to t - kj + 0 
t t t 


(4.19) 


Now from Eq. 4.18 the following expression can be derived: 


^ (cos a 4- j sin a) (cos 8 - j sin 8) = 

~ { cosa cos8 + sina sin8 + j (sina cos0 - sin8 cosa) } = 


r » - P 2 C 2 
r w + P 2 C 2 


= real 


(4.20) 
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This means that the imaginary part in the left-hand side of Eq. 4.20 
has to be zero: 

tan a - tan B (4.21) 


or: 


w. t + k. fc + 0. = a) t - £ + 0 

tr tr tr t t t 


(. 4 . 22 ) 


When the assumption is made: 


(0 — 0) — CO- 

t tr 2 


then Eq . 4.22 becomes : 


0 ^ - 0 „ 

t tr 


2 0 J 2 S, 


(4.23) 


From this expression follows that the phase difference between the 
transmitted wave and the transmitted-reflected wave is determined 
by the frequency of the sound wave and the cavity depth. This occurs 
when the speed of sound of air is assumed to be a constant. 

At z - 0 the impedance of the incident and reflected wave is 
identical to the sum of the normal specific acoustic impedance of 
the plate and the impedance of the transmitted and the transmitted- 
reflected wave: 


p. + p 

“i *r 


u. + u 

i 


r 


_ P-. + P 

z + — 

P u. u 


tr 


tr 


(4.24) 


where z^ represents the panel impedance. This impedance will be 
written as follows: 
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(4.25-) 


The Eqs. 4.12 and 4.13 already give the particle velocity of the 
transmitted and the transmitted-reflected wave. The incident wave 
and the reflected wave are also assumed to be plane progressive 
waves. Then the particle velocity of the incident wave is: 


u. 


p i c i 


C4.26) 


while the particle velocity of the reflected wave may be written 
as follows: 


u 

r 


p i c i 


(.4.27) 


Substitution of Eqs. 4.12, 4.13, 4.25, 4.26 and 4.27 into Eq. 4.24 
results in the following expression; 


or: 


P ± + P r p t + p tr 

= ~ p l c l ‘ r p + P 2 C 2 

p i - p r P t - p tr 


C4.28) 


p. + p 

r i r r 


P- - P 
r i *r 


p . c. = r + jx + 

11 P J P 


1 + p tr 7 p t 

- - P 2 C 2 

1 ' p tr 7 p t 


(4.29) 


The ratio of the transmitted-reflected sound pressure to the trans- 
aiitted sound pressure at z = 0 is: 



M« tr t + e ) 

D e 

j(fu t + 0 ) 

C e 


(4.30) 
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With the assumption and the use of Eq. 4.16, Eq. 4.30 

becomes : 


P 


tr 


p t 


D e 
C e 



w 


w 


- P 2 C 2 

+ P,c 2 


-2jk A 
e z 


04.31) 


Substitution of Eq. 4.31 in Eq. 4.29 results in the following expres- 
sion: 


p i + p r 

Pi c i = r + jx + 

p . - p 

1 *r 


-2jk„A -2jk ? A 

r w Cl + e * )+p 2 c 2 (l-e ) 


r (1 - e 
w 


-2jk 2 & 


-2jk 2 £ P 2 C 2 
)+P 2 c 2 (l+e ) 

(4.32) 


or: 


p. +p 


p,c, = r + jx + 

— — 11. p J r, 


r w +j p 2 c 2 tan k 2 A 


p . - p 
r i r r 


*p ‘ p 0 c_ + j r tan k 0 A P 2 C 2 
z z w z 


(4.33) 


The right-hand side, RHS, of Eq. 4.33 may be written as follows: 


(1 + can 2 k,)l)(p ? c ,) 2 r 
™ - * P * ; 2 . 2 2 2 2 ' + 


X + 

p 


( p 2 c 2 > + r w tan k 2 a 

2 2 

{ (P 2 c 2^ “ r w * P 2 C 2 tan k 2* 


(p 2 C 2 ) ' 


2 2 ,. 
r tan k 0 A 
w 2 


(4.34) 


or: 


RHS = k + j X 


(4.35) 


Eq. 4.33 may now be written as: 


p . + p 

r 


p. - p 

l r r 


p^ = k + 3 X 


(4.36) 
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From Eq. 4.36 follows: 


p jc — P-. c 1 + i X 
_ r _ K 1 1 J 

F. ’ c + °i c i + J x 


C4.37) 


By taking the magnitude of both sides of Eq. 4.37 (in the same way 
as done in Appendix A) the following result is obtained: 


B 

A 


, % 2 . ,2 

(k - P^) + * 

(k + J 


1/2 


(4.38) 


The sound power reflection coefficient is : 


a 

r 



(tc - p .^) 2 + X 2 
( K + P^cp 2 + X 2 


(4.39) 


where k and X follow from Eq. 4.34. 

At z = 0 Eq. 4.37 may also be written as follows (w^ = to = co^) : 



k - + j X 

K + + j X 


(4.40) 


where 0^ represents the phase angle between the reflected wave and 
the incident wave. From Eq. 4.40 follows: 



K 2 - (p 1 c 1 ) 2 + X 2 + jBAp^ 
(k + p-^) 2 + ^ 2 


(4.. 41) 


With this expression the phase angle 0 can be derived: 


36 



arctan 


(4.42) 


6 


r 


( 2 x p i c i ) 

( K 2 - (P 1 Cj _) 2 + X 2 \ 


At z = 0 also the following boundary condition must be satisfied: 
the particle velocities normal to the surface of the plate are equal. 
This is true when the panel damping is zero. This boundary condition 
can be written as follows: 


u. + u = u + u (4.43) 

i r t tr 

Substitution of Eqs. 4.12, 4.13, 4.26 and 4.27 in Eq. 4.43 results 
in the following expression: 




P t " P tr 
P 2 C 2 


(4.44) 


Now the assumption will be made: 


i. 


03 - 0 ), 


0 ). 

tr 


03 


In that case Eq. 4.44 becomes: 


— ^ 1 C 1 — — 
A - B =* (C - D) 

P 2 C 2 


(4.45) 


or: 

C = P 2 C 2 (1 - B/A) 
A P 1 C 1 (1 - D/C) 


(4.46) 


The ratio B/A follows from Eq. 4.37, while D/ C comes from Eq. 4.31. 
Substitution of these expressions in Eq. 4.46 results in: 
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(A .47) 


C 

A 


P 2 C 2 

p l c l 


2p l c l 


K + P 1 c 1 + j A 


r w CL 


- e 


-2jk 2 & 


) + p ? c 2 (l.+ .e 


-2jk 2 il 


T w + P 2 C 2 


or: 

— 1 + j tan k^Z 

A . P 2 C 2 ^ r w + P 2 C 2"* (k + p c, + j A) (p„c 0 + jr tan k„J.) 

±1 z z w z 

(.4.48) 


The denominator of the right-hand side of Eq. 4.48 may be written 
as follows: 


(kp 2 c 2 + p i c i p 2 C 2 “ Ar w tan k 2 ^ + 

j(xr ^ tan k^Z + p-^c^ tan k^Z + P 2 C 2 ^ (4.49) 


Now. the magnitude of Eq. 4.48 may be taken. This’ leads to the 
following result: - 


a“ p 2 c 2 (r w +p 2 c 2 ) - 


J= 

7 (KP^ + P^J 


1 -J* tan k 2 fc 


2 2 
p 2 c 2~* r w tsn + C<r^tan k^+p^c^) 


(4.50) 


I " 0 2 c 2- fr w +p 2 <: 2 ) ' 


1 + tan k 2 & 


(4.51) 


{ (k + P^ c ^) + A 2 )'{r w 2 tan 2 k 2 5. + (p 2 c 2 ) -2 } 


With above expressions the noise reduction through a cavity-backed 
flexible plate can be determined. Noise reduction is the difference 
in sound pressure level at the source side of the panel and the. 
receiver side, or: 
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NR = SPL 1 - SPL 2 


(4.52) 


where SPL^ is the sound pressure level at the source side and SPL, 
the sound pressure level at the receiver side of the panel, inside 
the cavity. Sound pressure level is defined as: 


SPL = 10 



(4.53) 


where P e represents the measured effective pressure (rms) of the 

sound wave and P c is the reference effective pressure. Here the 
ret 

rms sound pressure at the source side is P^, while the rms sound 
pressure at the receiver side is P 2 - Now Eq. 4.52 may be written 
in the following manner: 


NR = 10 log 



(4.54) 


The acoustic pressure p^ at a point at the source side of the 
panel is given by the sum of the real part of the incident sound 
wave p. and the real part of the reflected sound wave p . There— 

• L 

fore: 

p =* A cos (w.t - k. z ) + B cos (w t + k Z-, + 0 ) (4.55) 

1 i x 1 r r x r 

where 9^ indicates the phase difference between the incident wave 
and the reflected wave. In the same manner as done in Section 3.1, 
the rms pressure can be derived. 

When the frequency of the incident wave is not the same as the 
frequency of the reflected wave (uk t co^.) , then the rms pressure at 
the source side may be written as follows: 


39 



In the case that the frequencies of the incident wave and the reflected 
wave are identical (ok — , the rms pressure is r 


p i = 


f A 2 + + 2’ AB cos (8 + 2 z^) 


(4-57) 


This expression indicates clearly the importance of the microphone 
position. 

At the receiver side of the panel (inside the cavity) the 
acoustic pressure p^ at a fixed point tne sum the rea ^ P art 

of the transmitted sound pressure and the transmitted-reflected 
sound pressure. Therefore: 

P 2 = C cos (u) t t - k t z 2 + 9 t ) + D cos (w tr t + k tr z 2 + 9 tr ) 

(4.58) 

where 8^ indicates the difference in phase between the transmitted 
wave and the incident wave, while 0^ indicates the phase difference 
between the transmitted-reflected wave and the incident wave. In 
the case that the frequency of the transmitted sound wave is not 
identical to the frequency of the transmitted-reflected sound wave, 
no standing waves will occur inside the cavity and the rms sound 
pressure is not a function of the microphone position as shown by 
the following expression: 


■V 


2 P 

C + D 


(to t * 


(0 ) 
tr' 


(4.59) 
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When the frequencies of the transmitted sound wave and the transmitted- 
reflected sound wave are identical, then the rms sound pressure is: 

(4.60) 

Substitution of Eqs. 4.57 and 4.60 in Eq. 4.54 gives the following 
result : 

A 2 + B 2 + 2 A B cos (0 + 2 k z ) 

NR = 10 log ; (4.61) 

C + D + 2 C D cos (0 t - 0 - 2 k 2 z 2 ) 

B 2 B 

2 1 + —+2 j cos (0^ + 2 k x Zj) 

NR = 10 log ■ ■■ • y (4.62) 

C 1+S 2 + 2 § cos (e t‘ S tr" 2 k 2 z 2> 

C i 

where B/A follows from Eq. 4.38, 0 r from Eq. 4.42, D/C from Eq. 4.17, 

0 - 0 from Eq. 4.23 and C/A from Eq. 4.51. 

t tr n 

Eq. 4.62 is valid when the frequency of the incident sound wave 
and the reflected sound wave are identical and the frequency of the 
transmitted wave is identical to the frequency of the transmitted- 
reflected wave. However, to derive the expression for the ratio 
C/A, the assumption is made that the frequency of the sound waves 
at the source side is the same as the frequency at the receiver side 
or: 

to. - to = w = ok -to 
r r t tr 

In the following sections expressions will be derived for the 
impedance of the panel and the backing wall. 
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4.2 Panel Impedance of a Cavity-Backed Flexible Panel with Simply 

Supported Boundary Conditions 

The derivation to obtain the normal specific acoustic impedance 
of a cavity-backed flexible panel is basically identical to the one 
in Section 3.2, the impedance of a free flexible panel. In this 
section, too, the structural damping of the plate will be ignored. 
According to Ref. 14 structural damping usually has little or no 
effect on the natural frequencies and on the steady-state amplitude 
of the plate. 

The governing differential equation of motion of an isotropic 
flat plate is: 

D V 2 V 2 w (x, y, t) + m 8 v (x, y, t) . = (X} y> t) 

(4.63) 

This equation is based on the small deflection plate theory. 

In the following sub-sections the free vibration and the forced 
vibration of a simply supported cavity-backed flexible plate will 
be discussed. 

4.2.1 Free Vibration of Simply Supported Plate 

When the plate is vibrating freely, the external force, p z , 
is omitted and the differential, equation of motion is: 

D V 2 V 2 w (x, y, t) + m - — - - = 0 (4.64) 

S 3t 

This expression is identical to Eq. 3.43. Therefore, the 
solution will also be identical. Eq. 3.54 gives the natural fre- 
quencies of a freely vibrating plate, namely: 
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(4.65) 



where m-1, 2, 3, . . . ; n = 1, 2,3 . . . ; D represents the 
bending stiffness and the surface mass of the plate. 

4.2.2 Forced Vibration of a Cavity-Backed Simply Supported Plate 
The governing differential equation of motion of the plate is 
given by Eq. 4.63. The assumption will be made that the plate will 
be hit by plane waves only. In that case the forcing function is 
not a function of the coordinates x and y any more, but only a 
function of time. When discussing the forcing function of the 
free panel, it was concluded that the forcing pressure was the sum 
of the incident, reflected and the transmitted sound pressure. 

However, a cavity-backed panel is also excited by the transmitted- 
reflected wave or: 

P z (t) = (p ± + P r ) - (p fc + P tr ) (4.66) 

The differential equation of motion is linear. This means that 
Eq. 4.63 can be solved separately for p^, p^, p^ and P tr * In that 
case the total solution is the sum of the partial solutions. The 
total plate deflection can be divided in the transient deflection 
and the steady-state deflection w p . Here only the steady-state 
deflection will be considered. The total steady-state deflection 
of the plate will be: 

w_ ” w. + w 4 - 4 - (4.67) 

P i r t tr 
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In the same manner as done in Section 3*2.2 for the free panel, 
the steady-state deflection of a cavity-backed panel can be written 
as: 
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(4.68) 

(4.69) 

(4.70) 

(4.71) 

(4.72) 


With above expressions the normal specific impedance of the panel 
can be obtained. 

4.2.3 Normal Specific Acoustic Impedance of a Cavity-Backed Panel 

The normal specific acoustic impedance of a plate is defined as 
the ratio of the driving acoustic pressure to the resulting velocity 
of the plate: 
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(A* 73) 


z 

P 


p. + p - p_ 
x r t 


“ P 


tr 


w r 


The panel velocity can be obtained from Eq. 4.68 and may be written 
as follows: 


W P = ^“i w i * j“ r W r + ^“t W t + ^“tr w tr 


(4.74) 


where w^, w r> and w are given in Eq. 4.68. The panel impedance- 
may be written as: 


z 

P 


jui.t ;K<d t + 6 ) j(u t + 9 ) jCw t+0^) 

Ae +Be -Ce - P e 


j (ok w ± + eo r w r + m t w fc + m tr w tr ) 


(4.75) 


Substitution of Eqs. 4.68, 4.69, 4.70, 4.71 and 4.72 in Eq. 4.75 results 
in a very complicated expression. However, the assumption: 


(0. = 0) ~(J0 = 0) =0) 

x r t tr 


simplifies this expression and makes it more -useful. Because of this 


assumption the panel impedance is: 


; := ± (4.76) 

■n CO 00 

y i 16 „ sin(imrx/a) sin(mry/b) 

— -~ 2 ^ x l 5 2 

m it m~l n=l inn (to ^ - w ) 
s mn 


This expression is identical to the expression of the panel impedance 


of the free panel (see Eq. 3.86). 


The panel resistance r and reactance x are: 

P .P 
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r = 0 
P 


(4.77) 


x 

P 


_ r sin (mTrx/a) sin (n.ry/b) 

4 4 2 2 

m=l n=l m n (to - to ) 
mn 


These two expressions are identical to Eq. 3.88. 
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4.3 Normal Specific Acoustic Impedance of the Backing Wall 


In Ref. 15 (page 448) is mentioned that the simplest and most 
direct procedure to introduce cavity damping is to consider the 
absorbing characteristics of the wall as being determined by its 
normal specific acoustic impedance, z 9 the ratio of pressure to 
normal velocity at the surface of the wall. Ref. 15 expresses the 
normal specific acoustic impedance of a material relative to air 
by the expression: 

\ " <r n + * x n> P 2 C 2 ( 


Many of the materials, however, used for wall surfaces are charac- 
terized by r n » and r R » 1. The imaginary part of the wall 
impedance may be omitted, in most cases. 

In the model used in Ref. 9, the absorption characteristics are 
modeled by an impedance z . For representative acoustic materials: 


z 

w 


p 2 c 2 


= 1-10 


(4.79) 


In Ref. 9 also experimental results are presented for the damping 
ratio of the cavity as function of the frequency. The cavity wall 
absorption causes the damping. The curve obtained for the damping 
ratio of the cavity suggests that the cavity damping ratio at a 
constant cavity depth varies as the inverse of the frequency squared. 
This result, as mentioned by Dowell, is inconsistent with the theory. 
No explanation can be given for this damping result, other than 
this behavior indicates a more complex damping mechanism than that 
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of the panel. Also in this reference the wall reactance is assumed 
to be much smaller than the wall resistance. 

In Ref. 10 the cavity has an absorptive wall, too. The impedance 
of the absorptive wall is assumed to be f>2 C ‘2 at **igh frequencies. 

At low frequencies, however, the absorption material does not have 
any effect and the normal specific impedance of the wall will approach 
infinity as the frequency tends to zero. Therefore, a real relation- 
ship for is assumed that varies from infinity at zero frequency 
to p 2 c 2 at frequencies: 


w 


P 2 C 2 * 1 + ( 


100 

2(0 


> > 


(4.80) 


where represents the first cavity depth mode: 


ir c„ 


to 


100 - 


where % represents the cavity length (depth) . 


(4.81) 
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CHAPTER- 5 


RESULTS AND COMPARISONS 


In Chapter 3 the transmission of sound through a free panel 
has been discussed, while in Chapter 4 the transmission of sound 
through a cavity-backed panel has been examined. As can be expected, 
these two cases do not stand on their own. The noise reduction through 
a free panel can be seen as a special case of the noise reduction 
through a cavity-backed panel. The cavity behind a free panel 
can be considered to have an infinite size, or the backing wall 
has the same impedance as air: 

r w = P 2 C 2 C5 ' 1} 


where p 2 c 2 called the characteristic impedance (resistance) of 
the air at the receiver side. Substitution of Eq, 5.1 in Eq. 4.17 
results in: 

(5.2) 

or the pressure amplitude of the transraitted-reflected wave is 
zero, which means no reflection from the backing wall. Substitution 
of Eq. 5.1 into Eq. 4.39 gives the following result for the sound 
power reflection coefficient: 


a 

r 



^ r p + P 2 C 2 ~ p l c l^ + X p 
(r p + p 2 c 2 + p lCl ) 2 + x p 2 


(5.3) 


which expression is identical to Eq. 3.14. The same happens when 
Eq. 5.1 is used in Eq. 4.42, the expression for the phase angle: 
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arctan { 


} 


(5.4) 


6 = 
r 


2 x p-,c, 
p 1 1 


(r p + P 2 C 2 )2_ (P 1 C 1 )2 + X p 2 


This equation is the same as the one derived for a free panel, 

Eq. 3.17. As a last check Eq. 5.1 can be substituted into Eq. 4.51. 
This results in the following expression for the sound power trans- 
mission coefficient: 


a 


t 


c 2 ^cp 2 

a2 < r p + P 2 C 2 + p l c l )2 + x p 2 


(5.5) 


which expression is identical to Eq. 3.23, the sound power transmis- 
sion coefficient of a free panel. 

With Eq. 3.24 the transmission loss of a free panel can be 
computed. In Figure 5.1 the transmission loss for an 0.635 mm. (0.025 

in.) thick aluminum panel is shown. The size of the panel is 0.4572 
o 2 

x 0.4572 m.“ (18 x 18 in. ). The number of panel modes considered 
is seven. The material values given in Kinsler and Frey (Ref. 15) 
are used. The fundamental resonance frequency of the simply supported 
panel will be at 14.98 Hz. At that frequency the transmission loss 
of the panel will be zero. This will also happen at the higher 
odd, odd panel modes. At high frequencies, in the mass-controlled 
region, the transmission loss increases by 6 db for a doubling of 
frequency, as predicted by the "mass-law" (Ref. 13): 


TL = 10 log [1 + (- 


urn 
s 

2p 2 c 2 


) ] 


(5.6) 


According to this expression an 0*635 ram. thick aluminum panel will 
have a transmission loss, TL, of 28.3 db at a frequency of 2000 Hz. 
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The value computed with Eq. 3-24 compares well with the above value, 
as shown in Figure 5.1. Not shown in Figure 5.1 but demonstrated 
in Table 5.1 is the behavior of a free 1.016 mm. (0.040 in. ) thick 
aluminum panel at low frequency. 

Table 5.1: Sound Transmission through an 1.016 mm. Thick 

Aluminum Panel 


FREC (HZ) 


"L. (DB) 

RW. (RAILS) 

i °r 

2 

1 7 . SO 


- 2 .Cf 

rT.st 

11.75 

1 * 15.03 

r pf! 

1 t> on 

^ Cl 

n 1 c rw 

*.C0 

11.57 


h15.c: 

^ . on 

C . '1 

■nm 



ORIGINAL PAGE IS 
OF POOR QUALITY 

Low frequency means that the frequency is far below the fundamental 
resonance frequency of the panel. At low frequencies, in the stiffness- 
controlled region, the transmission loss will increase- by 6 db for 
each halving of frequency. This behavior, too, is well predicted 
by the method of Chapter 3. 

The noise reduction prediction which can be obtained with the 
expressions of Chapter 4 are compared with the results of the pre- 
diction method of Guy and Bhattacharya (Bef . 8). This method is 
briefly discussed in Chapter 2, while the expressions which are 
derived in Ref. 8 are recalled in Appendix D. In Appendix E a 
Fortran IV time-sharing computer program of the prediction model 
of Guy and Bhattacharya is given. This model, however , is only valid 
for perfect reflecting cavities; thus, cavity damping is not considered. 
To make a comparison between the two models possible; 
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(5.7) 


r = 1,000,000 kg/m^.sec 
w 

has been entered in the model -of Chapter 4. In that case the backing 
wall can be considered to be- a perfect reflector. In Appendix F a 
Fortran- IV time-sharing program of the prediction model discussed 
in Chapters 3 and 4 is given. By using both computer programs. 

Figure 5.2 has been generated. Figure 5.2 shows noise reduction 
curves of the following model: 

Panel: thickness Q.635 mm. (0.025 in.) 

material aluminum - 

length along x-axis 0.4572 m. (18 in.) 

length along y-axis 0.4572 m. 

Cavity: length along x-axis 0.4572 m. 

length along y-axis 0.4572 m. 

length along z-axis 2.7 m. (106.3 in.) 

The source and the receiver microphones 'both are situated at the 
axial-axis of the cavity. The source microphone position is 0.01 m. 
in front of the panel, while the receiver microphone is placed 
0.216 m. behind .the panel, inside the cavity. Also in this case 
the material values given in Ref. 15 are used. For both models 
7 unsymmetrical (odd, odd) panel modes are considered, and for the 
model of Guy and Bhattacharya 2 cavity modes are taken into account. 

Especially above the frequency of 100 Hz the agreement between 
the results of both models is very good. Below the frequency of 
100 Hz, the results of the model discussed in this- report show fair 
agreement with the results of the model of Guy and Bhattacharya (Ref. 8). 
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In Ref. 8 a' comparison is shown between experimental results, results 
of the model of Ref . 8 and results -of the prediction, method of Ref% 5. 
These results show good agreement with each other. This means that 
the behavior of a panel backed by a perfect reflecting cavity is 
fairly well predicted by the model of Chapter 4. 

In Figure 5.3 an experimental noise reduction curve is shown 
of a lead vinyl panel. The flexural rigidity of lead vinyl is very 
low, which means that the fundamental resonance frequency is very 
low, or: 


w 1:L “ 0 (5.8) 

In Figure 5.3, however, the noise reduction curve shows a major dip- 
at 45 Hz. This frequency indicates the first axial cavity mode. 
Substitution of this frequency in the following expression: 



results in a value for the "effective" cavity length. The "effective" 
-cavity length appears to be & - 3. 8 m. , which value is quite different 
from the actual cavity length l = 2.7 m. Around the frequency of 
700 Hz, severe dips and peaks start showing up in the noise reduction 
curve. This indicates that the cavity modes in x- and y-directions 
(see Figure 4.1) start showing up. The frequency of 750 Hz indicates 
an 0, 2, 0 cavity mode. This is the first cavity mode which will 
show up because of the position of the receiver microphone. 

In Figure 5 .,4 an experimental noise reduction curve is compared 

( 'i * 

with results of the prediction model of Chapter 4. The panel/cavity 
model used is almost the' same' as the one mentioned before in this 
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section. The only differences are that effective cavity length is 
used instead of the real length and that the impedance of the backing 
wall is a function of the frequency: 


r 

w 


co 


P 2 C 2 


< 1 + ( 


100 

2o) 


2 

) ) 


( 5 . 10 ) 


The agreement between the experimental and the predicted noise reduc- 
tion curve is poor. One cause for this poor agreement in the un- 
known behavior of the absorbing materials in the receiving chamber. 

An attempt has been made to model this behavior with Eq. 5.10. However, 
the results generated with this expression, but also with other expres- 
sions, compare poorly with the experimental results. A second cause 
is that the panel resonance frequencies show up too strong. Adding 
panel damping to the model will weaken the influence of the panel 
resonances. A third cause is that the boundary conditions for the 
model are simply supported, while the boundary condi tons for the 
experimental results are somewhere between simply supported and 
clamped. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


In this report a method has been discussed to predict the noise 
reduction through a cavity-backed flexible plate. The transmission 
of sound through a panel backed by a finite space is completely dif- 
ferent" from the transmission of sound through a panel backed by an 
infinite space* Figures 5.1 and 5.2 provide a good indication of 
this effect. 

A comparison of the results obtained with the prediction method 
of Guy and Bhattacharya (Ref. 8) and the results obtained with the 
method of this report show good agreement. Experimental results ob- 
tained with the test facility of Ref. 1, and results of the prediction 
method of Chapter 4 show poor resemblance. This, however, must be 
attributed partly to the experimental results. 

In the first place the boundary conditions of the test panels 
in the test facility are not known exactly. The edge conditions 
are somewhere between clamped and simply supported. Secondly, the 
behavior of the absorbing materials in the receiving chamber is not 
known. Thirdly, the cavity is not sealed airtight. A small air leak 
can change the panel behavior considerably. 

In the prediction method no panel damping has been included. 
Although the internal damping in thin aluminum and steel plates is 
not much, it has its effect on the size of the peaks and dips in the 
noise reduction curve. Damping decreases the size of the peaks and , 
dips. According to Figure 5.4 the influence of the panel resonances 
is too great. The use of the panel impedance equation derived by 
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Barton (Ref- 10) will omit the effect of the panel resonances in 
the noise reduction curve. Better results may be expected then. 

The prediction method derived in Chapters 3 and 4, however, 
is very useful in the investigation of 1) the effect of a change 
in the angle of incidence of the sound waves on the transmission 
of sound through a panel, and 2) the effect of a difference in air 
■pressure across a cavity-backed panel. In the first case, Kinsler 
and Frey (Ref. 15) indicate how to incorporate the angle of inci- 
dence of the sound waves into the expressions of Chapters 3 and 4. 

In the latter case, the normal specific acoustic impedance of the 
flexible panel also becomes a function of the pressure differential 
across the panel. The forcing function in the differential equation 
of motion of the plate will consist of a sound pressure term and an 
air pressure term. In Ref. 16 the effect of internal pressurization 
is discussed in the case of a curved plate. 
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APPENDIX A 


DERIVATION OF EQ. 3.14 


The following expression is known at z = 0: 


p r (r p * P 2 C 2 ~ °l' : l > * J *p 

F. (l p + P 2 C 2 + p l c l> + * x p 


This expression may also be written as follows: 


j 

B e 


(u t + 6 ) 
r r 


A e 


(r P + P 2 C 2 ~ p l c l ) + j X p 
(r p + p 2 c 2 + p lC]L ) + j x p 


(A.l) 


(A. 2) 


or: 

B (cos a 4- j sin a) _ a + jb 
A (cos g + j sin 6) c + jb. 


(A. 3) 


where: a = to t + 0 

r r 

g “ OKt 

8 - r p + P 2 C 2 ' p l c l 

b = x 
P 

c = r p + P 2 C 2 + p l c l 


(A. 4) 


From Eq. A. 3 the following equation can be derived: 
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A- 


(a 4- jb) (c - jb) (cos gjjyj sin ft) (cos a - j sin a) 


2 

c 


+ b 


2 


(A. 5) 


or: 


13 {ac-t-b + .jb(c - a)Hcos(a- ft) - i sin(a-g)} 

A 2 .I 2 

c + b 


(A. 6) 
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i>|w 


which results in: 


2 2 
(ac + b )cos(a-g) 4- b (c-a) sin(a-B) 4* j{b (c-a)cos (a-B)-(ac + h )sin(a-8)} 

2 2 
c + D 

(A. 7) 

Taking the magnitude of the right-hand side^ of Eq. A- 7 leads to the 
following result : 


B 

A 


2 -,2 
c 4 b 


(ac Jr b^) + (be - ab)^ 


or: 

B */Z±Z 

A " M c 2 + b 2 


Substitution of Eqs. A. 4 in Eq. A. 9 results in: 


B 
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+ p 2 c 2 - n^) 2 + Xp 2 

. . > 2,2 
+ p 2 c 2 + p^) + x p 


(A. 8) 


(A. 9) 


(A. 10) 
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.APPENDIX B 


TIME AVERAGE OF COSINE PRODUCT: 

Different Frequencies 

In this appendix it, will be shown that -the time average of 
the function: 

f(t) = 2AB cos (w^t - k^ z^) cos (<i> r t + k^ z^ + 0^) (B.l) 

is zero, if w^. The time average of a function is given as: 


+T 

<f> = lim x— =• / f(t) dt 

rrt t ^ ^ m 


(P.2) 


T «o -T 


The integral: 


+T 


/ cos (a), t-k.z, )cos (w t + k z, + 0 )dt = 
1 x 1 r r 1 r ' 


+T 

/ {cos w.t cos k.z, -i- sin w.t sin k.z,} 

_ J T i i 1 i i 1 

{cos w t cos (k z, + 0 )-sinw t sin(kz, +9 )} dt = 
r rlr' r r 1 r 


4-T 


(1) cos k.z, cos(k z, + 0 ) / cos w.t cos w t dt + 

v ' il V rlr /J T i r 


+T 


(2) -cos k.z, sin(k z, + 0 ) / cos w.t sin w t dt + 

il r 1 r' J T x r 


+T 


(3) sin k.z, cos(k z, +0 ) / sin w.t cos w t dt + 
xl r 1 r *_ i r 


-T 


+T 


(4) -sin k.z, sin(k z, + 0 ) / sin w.t sin w t dt 

xl r 1 r _ T i r 


(B.3) 
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The solution of integral (1) is : 


+T 

j COS 01. t COS 01 t dt = 

J i r 

-T 


sin(w. +o) )t sin(oi.-o) ) t 
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2(^ + 0^) 2 (ok - oj^). 


+T 
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(B-4) 


The solution of integral (2) is: 


+T 


-T 


cos oi.t sin oi t dt = 
x r 


cos (uk- oi r ) t cos (ok + 0 )^) t 


-I+T 


2(o).-o) ) 
x r 


2(oi. + oi ) 
x r 


= 0 (B.5) 


J_ T 


In that case the solution of integral (3) is: 


+T 

/ sin OKt cos to r t dt = ‘0 (B.6) 


The solution of integral (4) is: 


4-T 


/ sin oi. t sin oi t dt = 
-T 
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2 (a) , 4* to ) 
x r 


-i+T 


U T 
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(B.7) 


The solution of Eq. B*3 becomes: 


sin(a3.-hoj )T sin (ok- &> r )T 

cos k.z- cos(k z- + 0 ){ 7 7 7 r 
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(B.8) 
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According to Eq. B.2, Eq. B.8 has to be divided by the factor 2T. 

Then the limit for T -*■ 00 has to be taken: 

<f> = lim •^jr*(Eq. B.8) = 0 (B.9) 

T ->• «> ZL 
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Identical Frequencies 
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APPENDIX C 


TIME AVERAGE OF COSINE PRODUCT: 

Identical Frequencies 

The time average of the function: 

f(t) = 2AB cos (dKt --k^z^) cos (a^t + k^z^ + 9^.) 

is not equal to zero when: 

0). = 0) = (I), 

i r 1 
Then Eq. C.l becomes: 

f (t) = 2AB cos (w^t - k^z^) cos (co^t + k^z^ + 0^) 

The time average of f(t) is defined as: 

1 +T 

< f > - lim / f (t) dt 

T -> co -T 

The integral : 

+T 

/ cos(oj^t - k^z^)cos(a)^t + k^z^ + 8 r )dt = 

+T 

/ {cos d^t cos k^z^ + sin u^t sin k^z^} 

{cos &>^t cos (k^Zj + © r ) - sin w^t sin(k^z^ 

+T 2 

(1) cos k,z- cos(k. z. + 6 ) / cos w, t dt + 

11 11 r _t 1 

+T sin 2u), t 

(2) -cos k^z^ sin(k^z^ + 0^) / 2 dt + 


(C.l) 


(C.2) 


(C.3) 


e r )}dt » 
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(3) 


dt + 


-FT sin 2co^t 

sin k- z-. cos(k n z. + 0 ) / 

11 11 r Jx 2 


-FT 


(4) -sin- Lz- sin(k.z- 4* 0 ) / sin a)., t dt 

11 11 r -x 1 


The solution of integral (1) is: 



The solution of integral (2) and (3) is: 


-FI sin 2a)- 1 
f 1 


cos 2u)^t 


+T 


-T 


= 0 


The solution of integral (4) is: 


+T 


-T 


2 

sin co^t dt — 


sin 2w, t 
t i 

2 " 4(0, 


-i+T 


= T - 


sin 2w^T 

2 ( 0 , 


T 


The total solution of Eq. C.4 is: 


sin 2(o^T 

cos + 6 r ){T + — ) + 

sin 2(0- T 

- sin k^ sin(k lZl + 0 r >{T - — > 


The time average: 


<f > 


lim ~-*(Eq, C. 8) 
T ->• <*> 


(C.4) 


(C.5) 


(C.6) 


(C.7) 


(C..8) 


(C.9) 
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< f > = AB{cos k^z^ cos (k^z^ + 8^) - sin k^z^ sin(k^z^ + 0^.)} 



(C.10) 

AB cosCk-Z, + k. z +0 ) 
11 11 r 

(C.ll) 

AB cos (2 k. z. + 6 ) 
11 r' 

(C.12) 
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APPENDIX D 

NOISE REDUCTION PREDICTION METHOD 
BY GUY' AND BHATTACHARYA 
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APPENDIX D 


NOISE REDUCTION PREDICTION METHOD 
BY GUY AND BHATTACHARYA 

The model considered by Guy and Bhattacharya is shown in Figure 
D.l and consists of an acoustically hard walled rectangular cavity 
having at one face (x = 0) a flexible vibrating plate. 


y 



X- 


Figure D.l: Sketch of Panel and Cavity Geometry 


The following expression is derived for the acoustic pressure within , 
the cavity: 


P 7 (X, y, z. t)= Z K.K'. 

q=l 

r=l 


m=0 


n-0 


B [Z (0)cav]cosh[ (a-x)a(t) ] 
mn mn 

r, qr .ITITT . .H7T N jhJt 

lqr [2^(0)0317+ Z^ r pan]cosh(aa(t)] cos ^b ” COS c z ^ e 


(D.l) 
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where : 


K = 0.5 for m = 0 K = 1 for m > 0 


K'=0.5 for n = 0 


K*= 1 for n > 0 


The assumption that the pressure wave on the surface of incidence, 
P^, is a plane wave makes the following simplification possible: 


P x 16 


2 for q and r odd 


lqr 


qr ir 
0 otherwise 


(D.2) 


The coupling coefficient B may be written as follows: 

mn 

qr 


B 1 r 

mn 2 / 2 ^ v ^ v 

qr ir (m - q ) (n -r ) 


mn 2 — 2~2 2 — T~ [ cos (^)c°s(nnr)- [cos(rir)cos(mr)-l] 


The cavity impedance Z^(0)cav is given as follows: 


(D.3) 


^ . coth[aa(t)] 

a) > w: Z (0)cav « j p o) ttt 

mn — mn J o a(t) 


(D.4) 


where : 


a<t) 


2 r 

CO - ID 
mn 


(D.5) 


to < co: Z (0)cav = 

mn mn 


- j P o co 


cot[aa(t) ] 
a(t) 


(D.6) 


where : 


a(t) 


c 


co - CO 


mn 


CD. 7) 


According to Guy and Bhattacharya the following expression can 
be used for the natural frequencies of the cavity: 
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/ 2 2 
<k r m v . n v 
co = tic w ( 7~ ) 4- ( — ) 

ran 0 v b c 


(D.8) 


The expression for the panel impedance is: 


„ . ph , 2 2 N 

Z_, pan = j — (w - a) ) 
lqr r J (0 qr 


(D.9) 


where w is the panel natural frequency, which may be obtained from 
the following expression in the case of simply supported boundary 
conditions : 


<•„ - » 2 [ c ?■ ) + ( f > 1 


(D.10) 


where D represents the bending stiffness of the plate and p and h 
the density and the thickness of the plate. 

When the acoustic pressure within the cavity is measured at an 
arbitrary point on the axial centerline of the cavity, ^/2, c/2, t), 

then the term: 


cos ( y)* cos ( z) for m, n = 0, 1, 2, 3, . . . (D.ll) 

b c 


can be changed into: 


cos ( ~ ) cos ( ) for m, n = 0, 1, 2, 3, 


(D.12) 


m + n 

(-1) ^ for m, n = 0, 2, 4. 


(D.13) 


Substitution of the expressions mentioned above in Eq. D.l and by taking 


the magnitude of this expression, the noise reduction can be computed: 


77 



(D.14) 


1 

NR = 20 log 1 

2 

Appendix D gives a computer program of this me mod. . 


Symbol 

a 


b 

B 

mn 


c 

c 

o 

D 

h 

J 

K 


K T 


m 


n 


P 


1 


q 

r 

t 

x 


Table D * 1 : List of Symbols 

c 

Definition 

Room dimension in x~direction 
Room dimension in y-direction 
Coupling coefficient between room and 
panel modes 

Room dimension in z-direction 
Velocity of sound in air 
Modulus of rigidity 
Thickness of panel 

/T 

Constant 

Constant 

Room mode in y-direction 
Room mode in z-direction 
Acoustic pressure on the surface of 
incidence 

Acoustic pressure within the cavity 
Panel mode in y-direction 
Panel mode in z-direction 
Time 

Model coordinate 


Dimension 

m 

m 


m 

m/sec 

Nm 

m 


N/nT 

N/m^ 


sec 
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Table D.l:. List of Symbols (continued) 


Symbol 

Definition 

Dimension 

y 

Model coordinate 


z 

Model coordinate 


Z (O)cav 
mn 

Cavity impedance of an (m, n) cavity mode 
at the surface of the panel at the re- 



ceiver side 

kg/m“.sec 

Z oan 
qr F 

Panel impedance of the (q, r) panel mode 

kg/m~.sec 

P 

Density of panel material 

kg/m 3 

p o 

Density of air 

kg/m 3 

0) 

Forced angular frequency 

rad/sec 

0) 

mn 

Natural frequency for room mode (m, n) 

rad/sec 

0) 

Natural frequency for panel mode (q, r) 

rad/sec 
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APPENDIX E 


TIME- SHARING PROGRAM OF METHOD OF 
APPENDIX D 

In this appendix a Fortran IV time-sharing computer program 
will be given, which computes the noise reduction through a cavity- 
backed flexible panel as function of the frequency. The method 
used in this program is discussed in Appendix D, and the definition 
of the variables can be found in Table D.l. Table E.l provides a 
listing of the variables used in this program, while Figure E.l 
shows the flowchart. Figures E.2, E.3 and E.4 give the listing 
of the program. 


Symbol 

A 

ALPHA 

B 

BMNQR 

BNR 

C 

CA 

G0NST1 

C0NST2 

C0NST3 

COSHD 


Table E.l: Variable Names in Program 


Definition 

Room dimension a in x-direction 

Eqs. D.5 and D.7 

Room- dimension b in y-direction 

Coupling coefficient B 

rnn 

qr 

Noise reduction NR 

Room dimension c in z-direction 

Velocity of sound c q in air 

Constant K 

Constant K T 

16/ qrir 2 

Cosh[aa(t) 3 


Dimension 

m 

-1 

m 

m 


db 

m 

m/sec 
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Table E.l: Variable Names in Program (continued) 


Symbol 

Definition 

Dimension 

COSHN 

Cosh[(a - x)a(t)] 


D 

Bending stiffness D 

Nm 

DUM1 

(a - x)ci(f) 


DUM2 

a*a(t) 


E 

Young’s modulus E 

N/m^ 

F 

Frequency f 

Hz 

FACTOR 

Eq* D*1 without summation sign 


H 

Plate thickness h 

m 

ICMAX 

Number of cavity modes 


IPMAX 

Number of panel modes 


IHIGH 

Maximum frequency 

Hz 

ILOW 

Minimum frequency 

Hz 

ISTEP 

Increase in frequency 

Hz 

OMG 

Angular frequency to 

rad / sec 

OMGMN 

Natural frequency to for room mode 

rad/ sec 

OMGQR 

Natural frequency w for panel mode 

rad/sec 

PI 

ir 


Q 

Panel mode q in y-direction 


R 

Panel mode r in z-direction 


EM 

Room mode m in y-direction 


RN 

Room mode n in z-direction 


ROA 

Density of air p Q 

kg/m 3 

V 

Poisson 1 s ratio v 


X 

Receiver microphone position in x-direction 

m 
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Table E.l: Variable Names in Program (continued) 


Symbol 

Definition 

Dimension 

ZMNCAV 

Cavity impedance Z (O)cav 
mn 

2 

kg/m -sec 

ZQRPAN 

Panel Impedance Z^pan 

kg/m~*sec 
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Eq. D.3 
Eq. D.10 
Eq. D.8 



Eq. D.5 and D.7 
Eq. D.4 and D.6 
Eq. D.9 
Eq. D.l 


Figure E.l: Flowchart of Program 
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Figure E.2 : Listing of Program ORIGINAL 

OF POOR QUALITY 
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Figure E.3 : Listing of Program (Continued) 
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Figure E.4 : Listing of Program (Continued) 
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APPENDIX F 

TIME-SHARING PROGRAM OF METHOD 
OF CHAPTER 4 
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APPENDIX F 


TIME-SHARING PROGRAM OF METHOD 
OF CHAPTER 4 

In this appendix a Fortran IV time-sharing computer program 
will- be given to predict the transmission loss of a panel and the 
noise reduction through a panel. The panel can be either free or 
cavity backed (1=1 stands for cavity-backed; I £ 1 stands for 
a free panel). In the case of a free panel, some of the input 
data are not relevant and will not be used by the computer program. 
These data are: CL and Z2. Table F.l explains the different 

variables of the program. Figure F.l gives a flowchart of the 
program, and Figures F.2, F.3 and F.4 give a listing of the 
program. 


Table F.l: Variable Names in Program 


Symbol 

A 

ALABD 

ALPHAC 

ALPHAR 

ALPHAT 

B 

BNR 

CA1 

CA2 


Definition 

Room dimension a in x-direction 

Imaginary part A of Eq. 4.34 

D 2 /C^, see Eq. 4.17" 

? 2 

B /A , see Eq. 4.39 
C 2 /A 2 , see Eq. 4.51 
Room .dimension b in y-direction 
Noise reduction NR 

Velocity of sound in air at source side c^ 
Velocity of sound in air at receiver 
side c^ 


Dimension 

m 

kg/m 2 . sec 


m 

db 

m/ sec 
m/ sec 
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Table F.l: Variable Names in Program (continued) 


Symbol 

Definition 


Dimension 

CAPPA 

Real part k of Eq. 4*34 


kg/m 3 . sec 

CL 

Room dimension £ in z-direction 


m 

D 

Bending stiffness D 


Nm 

E 

Young’s modulus E 


N/m 2 

F 

Frequency f 


Hz 

H 

Plate thickness h 


m 

IHIGH 

Maximum frequency 


Hz 

ILOW 

Minimum frequency 


Hz 

ISTEP 

Step in frequency 


Hz 

MAX 

Maximum number of panel modes 



OMG 

Angular frequency 


rad/sec 

0MG1 

Fundamental cavity mode in z-direction 

rad/sec 

PI 

IT 



ROA1 

Density of air at source side 


kg/m 3 

R0A2 

Density of air at receiver side 


kg/m 3 

ROC1 

Characteristic impedance of air at 

source 



side 


kg/m 3 * sec 

ROC 2 

Characteristic impedance of air at 

receiver 



side p 2 c 2 


kg/m 3 . sec 

ROP 

Density of plate material p 


kg/m 3 

RP 

Panel resistance 


kg/m 3 - sec 

RW 

Wall resistance 


kg/m 3 .sec 

TM 

Tan 



THETAR 

Phase angle of reflected wave 0^ 


rad 
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Table F.l: Variable Names in Program (continued) 


Symbol Definition Dimension 

TL Transmission loss db 

X Receiver microphone position in x- 

direction m 

2 

XP Panel reactance kg/m .sec 

Y Receiver microphone position in y- 

direction m 

Z1 Source microphone position in z-direction m 

Z2 Receiver microphone position in z- 

direction ra 
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Figure F.2 : Listing o£ Program 


ORIGINAL PAGEK> 
OP POOR QUALITY 
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Figure F.3 : Listing of Program (Continued) 
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Figure F.4 : Listing of Program (Continued) 
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